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Abstract 

We present a class of theories of two dimensional gravity which admits homogeneous and 
isotropic solutions that are nonsingular and asymptotically approach a FRW matter domi- 
nated universe at late times. These models are generalizations of two dimensional dilaton 
gravity and both vacuum solutions and those including conformally coupled matter are 
investigated. In each case our construction leads to an inflationary stage driven by the 
gravitational sector. Our work comprises a simple example of the 'Nonsingular Universe' 
constructions of ref . [1] . 
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1 Introduction 



The singularity theorems of Penrose and Hawking^] pose serious physical and philosophical 
problems concerning the behaviour of physics in certain regions of spacetime. Taken at face 
value the theorems tell us that the manifolds described by General Relativity (GR) are 
generically geodesically incomplete - there are always geodesies that cannot be continued 
for arbitrary values of their affine parameter. This conclusion is highly unsatisfactory for a 
number of reasons. 

Firstly, the existence of singularities limits the predictive power of any physical theory 
for which the manifold is a background. Particles reaching the singularity cease to exist, 
carrying with them vital information about the particle's history. 

A second problem involves the fact that, as experience with 'physical' solutions to GR 
shows, the existence of a singularity is often accompanied by the unbounded increase of one 
of the curvature invariants of the theory. Examples of this include the infinite Weyl tensor 
at the centre of a Schwarzschild black hole and the 'big bang' singularity encountered in 
cosmological models. In fact it can be shown that if the singularity is reached on a timelike 
curve in a globally hyperbolic spacetime then the Riemann tensor Rfj, upa is unbounded ||. 

Finally, in connection with both of the above there is the problem of the quantum 
information loss paradox. This is associated with the evolution of pure quantum states into 
mixed states during the evaporation of black holes. 

It is a common belief that when we have a complete quantum theory of gravity the 
singularities predicted for GR will not occur. The appearance of singularities should be 
seen as a signal of the breakdown of GR and an indication of the need for a new theory. 

Since there is at present no theory of quantum gravity we shall attempt a more modest 
goal. We shall implement the 'Limiting Curvature Construction' (LCH)Q in a class of two 
dimensional theories of gravity. By doing so we shall obtain theories to which all solutions 
are singularity free and the only curvature invariant, the Ricci scalar R, is bounded. Thus, 
our key assumption is that those processes that may be possible for curing the 'sickness' of 
GR are representable at the classical level in a modified theory of gravity. 
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The models we will consider are generalizations of two-dimensional Dilaton Gravity 
which is obtained from the low energy limit of string theory || and is a theory of gravity 
coupled to a scalar field, the dilaton. We generalize these models by allowing an arbitrary 
coupling between the fields and a general potential for the scalar (which we still refer to as 
the dilaton). In a theory derived from string theory such a scalar field would be dynamical 
and the analysis of the model would be quite complicated. Here, however, we treat the 
scalar as a Lagrange multiplier. 

Nonsingular black hole solutions to this theory were studied in|| and so, in this paper, 
we shall concentrate on cosmological solutions to the theory. Cosmological solutions to 1+1 
dimensional dilaton gravity have been studied by other authors^] but our goals differ from 
theirs in that we wish to study generalizations of that theory with emphasis on the possible 
removal of spacetime singularities. 

The outline of this paper is as follows. In section 2 we review the model of general- 
ized dilaton gravity and obtain the field equations. In section 3 we explicitly search for 
homogeneous cosmological solutions and show that our construction yields theories with 
a 'natural' initial inflationary stage evolving into a FRW matter dominated universe. All 
manifolds described by this theory are singularity- free. In section 4 we investigate the effect 
of including a conformally coupled scalar field into the theory and show that this does not 
affect our ability to construct nonsingular cosmological theories. Finally in section 5 we 
conclude and speculate on the implications of our results. 

2 Two Dimensional Dilaton Gravity 

Here we shall describe the two dimensional models that we have investigated. Given that 
the model of dilaton gravity discussed in the previous section appears to be important for 
physics we shall study general two dimensional models describing a massless scalar field 
(which we shall still refer to as the dilaton) coupled to gravity. Our aim is to identify 
classes of theories to which all solutions are nonsingular. In this way we hope to pinpoint 
those features of the models which may be important for gravitational physics if nonsingular 
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behaviour is manifested at the classical level. We shall present the model and derive the 
general field equations in preparation for the implementation of our construction in the next 
section. 

Our starting point is the most general renormalizable lagrangian for gravity coupled to 
a scalar field in 1+1 spacetime dimensions. 

C = y/^[D(4>)R + G^XV^) 2 + F(0)] (1) 

with action 

S[g,4>}= f Cd 2 x (2) 

where D ((f)) is the coupling of the dilaton to gravity, G(4>) is just a conformal prefactor for 
the dynamical terms for the dilaton and H(<p) is a potential for the dilaton. 

We are free to perform a conformal or Brans-Dicke transformation on this Lagrangian[||] || 

via 



If we now require that 



g old = e 2r{<t>) g n 7 (3) 



~ = -m (4) 



then we may express the Lagrangian in the form 



C = ^{V{4>)+D{4>)R) (5) 

where the new potential is given by V{4>) = #(»e 2TW . This is the fundamental La grangian 
which we shall consider. In particular note that it is the metric appearing in this expression 
that governs the motion of test particles. Note that in principle we may eliminate 4> from 
this Lagrangian and thus obtain a higher derivative theory describing the gravitational field 



c new = jr^ (g) 



Now consider the equations of motion arising from the variation of the action with 
Lagrangian as in equation (5). Variations with respect to the scalar field yield 




(7) 



Now consider variations with respect to the metric 



5S = 



J d 2 x [(V(<P)+D(<p)R)^g a ^5g^ + 

{-D{<P)R aP 5g aP - D(<p)V 2 g aP 8g a P + D(<f>) V a Vp5g^)\ 



(8) 



Integrating by parts twice, ignoring surface terms and requiring 5S = gives 



D^Raf, + (V Q V^ - g a pV 2 )D{4>) + ^V(<f>)g afi + ^D(cf>)Rg a p = 



(9) 



Finally we note that in two dimensions R a p = g a pR/2 which puts our field equations in the 
form 



Equations (7) and (10) are our basic equations. In particular, (7) is the constraint equation 
due to the dilaton and it is through this that we shall force our theories to have limiting 
curvature and to be singularity free. In ref. 6 nonsingular static black hole solutions to 
this model were constructed. Here we seek to address the case when the spacetime is 
homogeneous and isotropic. To be exact we shall search for singularity free cosmological 
models. 

3 Cosmologies in Two Dimensions 

In this section we shall construct a vacuum model for which all homogeneous and isotropic 
solutions are nonsingular and approach a matter dominated compactified FRW spacetime 
at late times. Note that in particular these spacetimes exhibit no big-bang singularity. 



V(4)g a p = 2{V 2 g aP - V a V /3 )L»(0) 



(10) 
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We begin by making a redefinition of 4> so that D(4>) = We are perfectly free to do 
this and the choice is made for calculational simplicity. Since we are looking for spatially 
homogeneous solutions we shall use the ansatz 

g^ v = diag(-l,a 2 (t)) (11) 

for the metric. With these choices the field equations become 



v W -£*=o (i2) 



and 



where dot denotes differentiation with respect to time. Also our constraint equation (7) 
becomes 

= R (14) 

We may combine (12) and (13) and, for dcf)^ 1 /dt / 0, integrate to obtain the scale factor 
a{t) 

a(t) = - 7 A (15) 

where 7 > is an arbitrary integration constant. Substituting into (12) we have 

V((j)) = -7 _1 2a (16) 

Now, our construction consists of the following steps. Firstly we shall require that the scale 
factor approach that of a matter dominated FRW spacetime at late times. We shall then 
solve for the asymptotic behaviour of the potential in this region. We then demand that 
the curvature be bounded at early times. This yields a particular form for the scale factor 
which allows us also to solve for the asymptotic form of the potential in the early universe. 
Finally we shall choose a well-behaved interpolating potential and use this to solve for the 
general scale factor for our model. 

At late times, t — > 00, we shall impose a(t) = vt 2 ^ , where v > is a constant. Integra- 
tion of (15) yields 

m = \lt-^ (17) 



up to an integration constant. Prom (16) the late time form of our potential is given by 

vw = -| VV3 = _ a0 l/5 (18) 
37 

where a is a constant constructed from the other constants we have introduced. 

Now, at early times, t — > —00, we require the curvature, R, to be bounded, R = 20 1 
say. By (14) this implies 

-2B 2 

V(cp) = (19) 
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where (3 > is a constant. Then (13) becomes a harmonic oscillator equation for eft -1 giving 

<j){t) = JoT (20) 

y ' die! 31 + a 2 e~l 3t y ' 

with a\ and 02 arbitrary constants. Further, (15) can be rewritten as 

a(t)= 7 j t (cf > - 1 ) (21) 

So for various choices of a\ and 02 we can obtain an expanding de Sitter solution with 
Hubble constant H = 2f3 2 , (02 = 0), a contracting de Sitter solution (ai = 0) or a de Sitter 
bounce (ai = —02). 

Note that we have obtained an initial inflationary phase by bounding R, the only cur- 
vature invariant in 2 dimensions, and thus our solutions are nonsingular (both in the sense 
of avoiding infinite curvature and having a geodesically complete manifold) . It is an attrac- 
tive feature of this and other higher derivative theories of gravity that inflation (jUJ, which 
provides a solution to the homogeneity, horizon and flatness problems in cosmology, can be 
the result of the action of gravity alone and needs no matter fields to drive it. 

Now to find a model close to a compactified version of our real Universe, we shall search 
for a solution which has a de Sitter expansion at early times and evolves into a 'matter- 
dominated' FRW universe at late times. To this end we choose 

V Up) —>■ as (p — > 00 

9 

V{4>) -» -«(/> 1/5 as (p^O 
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A possible choice for the interpolating potential is 

V 'W = - l + a( 2ffl)-'^ <22> 
Explicit solution for the scale factor is not possible but it is only neccessary to show that 

our metric is well behaved as it interpolates between its asymptotic limits. From (13) we 
have (in terms of D = (j)^ 1 ) 

2D = (23) 

D 6 / 5 + b v ; 

where b = a/2/3 2 . This may be integrated to give 

D 2 = af(D) (24) 



where 



f(x) = |r 4 / 5 - ^ 6 2/3 ln(x 4/5_ 6 l/3 x 2/5 + 6 2/ 3) 

5 l2/3 ,2x 2 /5_6V3 
,b 1 ' 6 arctan( 7 =— — ) 



2V3 ' v^fe 1 / 3 

+ ^2/3 ln(x 2/5 +6 l/3 ) __^_ 6 2/3 arctan( ^3 ) (2g) 

The scale factor, a(t), may be expressed implicitly using (15) as 

a(D) = 7 D = 7 ^Jaf(D) (26) 

where we have chosen the positive square root for consistency. In Figure 1 we plot a(D) for 
a = 2/3 2 = 7 = 1. The analysis is not qualitatively sensitive to the value of these positive 
constants. 

It can be seen that a(D) is a real, positive, monotonically increasing function of D for 
D > 0. Further, when the universe emerges from an expanding de Sitter phase, -^D = 
D/D > for all D > and D is positive for all time because D — > + as t — > — oo and 
D(D) > for all D and so we can use the fact that D > V t to assert that the scale factor 
is also a monotonically increasing function of time. Our construction is therefore complete. 

In order to check the above results we also studied numerically the case of a contracting 
universe to ensure that the model made a smooth transition between the two asymptotic 



regions. Beginning with final conditions in the t — > oo, matter dominated region we evolved 
the equations backwards to the initial de Sitter solution at t — > — oo. Our numerical results 
confirmed the above analysis. 

Note that, although (p ^ oo as t — > — oo this is in fact a coordinate singularity. The 
coordinates t, x which we have chosen only cover half of de-Sitter space. If we consider 
de-Sitter space as a hyperboloid embedded in a flat three-dimensional space then we may 
define new coordinates which describe the whole of the space[ll|. In these coordinates it is 
clear that the point t = — oo is a coordinate singularity and that the space is geodesically 
complete since geodesies may be smoothly continued through the point t = —oo.. 

To summarize, we have shown that the class of theories with potentials satisfying the 
asymptotic conditions (18), (19) have all homogeneous isotropic solutions which are non- 
singular. Further, in bounding the curvature at early times our solutions have a 'natural' 
inflationary stage driven by the gravitational sector. Figure 2 shows a log-log plot of the 
scale factor demonstrating the interpolation between the behaviours of the asymptotic re- 
gions. 

We shall now turn our attention to modifications to the vacuum theory. 



4 Inclusion of Conformal Matter 

Now that we have found a class of theories for which all vacuum cosmological solutions 
are singularity free we wish to extend our analysis to consider the effect of mattter on our 
construction. 

In these two dimensional models the requirement of quantum consistency of the theory is 
equivalent to the conformal invariance of the Lagrangian. We plan to include matter using a 
scalar field and since the dynamical terms imply that a scalar field must transform trivially 
under a conformal transformation, we may not include a simple, polynomial potential in the 
Lagrangian. Thus, by far the most simple choice is to conformally couple a scalar field into 
the theory using just the dynamical terms. We shall represent this matter by a and shall 
demonstrate that the presence of this field does not prevent us from constructing singularity 
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free theories. We shall again restrict ourselves to search for theories with homogeneous, 
isotropic cosmological solutions. 
Our action is now 



S[g,<f>;a] = J {C g + C m ) d 2 x (27) 

where C g represents the gravitational sector of the theory given by the previous Lagrangian 
(26) and C m represents the Lagrangian for the matter described by the conformally coupled 
scalar field a. 



C m = X^d^a)^ (28) 

The equations of motion derived from this action are the constraint equation (14), the 
equation of motion for a 



V 2 <r = (29) 

and the modified field equations for the metric which, with the ansatz (11) and D(<j>) = (j)^ 1 , 
take the form 

and 

4--^(0) + ^ 2 = O (31) 



Since we seek solutions that depend only on time, we may solve (29) exactly to give 



& = - (32) 
a 

Substituting this into our field equations and combining we obtain 



h A2 " 2 =0 (33) 



b 2 a a 2 



and 
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Now, as before, let us search for solutions that asymptotically approach a matter-dominated 
FRW universe with scale factor given by a(t) = ut 2 ^ '. In this region it is straightforward 
to solve for the evolution of the dilaton and we obtain 



3 ^ 2 ° 2 ^2/3 _ ^5/3 



(35) 



2v 2 ' 5' 

where 5 is a constant of integration. This expression is valid for late times, t — > oo. 

For 5^0 the second term in eq.(35) dominates. Eq.(34) can then be solved for the 
potential at late times as 

V{(f>) -»• -^ 1/5 as t -> oo (<£->()) (36) 

where /i = (3<5/5) 1/5 (4(5/3) . 

At early times, t — > 0, we shall impose the condition that the curvature remain bounded, 
R = 2f3 2 , where 2(5 2 is a constant representing the limiting curvature. Equation (14) then 
implies 

-2B 2 

V(<f>) = —j— as t -> (37) 

9 

For the purposes of this paper we shall concentrate on the solution where the universe 
emerges from a de Sitter bounce at t = since this case turns out to be the easiest to treat 
analytically. Thus we write the scale factor at early times as 

a(t) = cosh(/3t) (38) 

where we have chosen a(0) = 1 without loss of generality. We can combine (33) and (34) 
into the following two equations, written in terms of D = 4>~ l 

K 2 

2D + V(D) + ^2=0 (39) 

2 

2D- + V(D) - ^ = (40) 
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where k 2 = \ 2 a 2 . Equation (39) with a(t) from (38), linearized about t = , admits the 
solution 



where c is an integration constant. Provided that 



D(t) = Trm - 2c cosh (pt) (41) 



(42) 



(40) is also satisfied by this solution to linear order in t. For cCl, D{t) tends to the small 
constant — c as t — > , so that the limit t — * corresponds to \4>\ ^> 1. To connect the two 
asymptotic regions of the de Sitter bounce and matter domination, we can then find a form 
of the interpolating potential similar to the one used in the vacuum model 

V ^> ~ " (1 + c0) 6 / 5 + ^(2/3 2 )~V 6/5 
It is an attractive feature of our theory that similar potentials can describe the nonsingular 

cosmological evolution of a universe with and without matter. We cannot solve for the time 

evolution analytically with this potential, so we evolve (14) and (39) numerically with a 

Runge-Kutta algorithm for a(t) and D(t), using the fact that in two dimensions with our 

ansatz for the metric, the Ricci scalar takes the simple form R = 2a/ a. In Figures 5,6 and 

7 we plot the scale factor a(t) , the Ricci scalar R(t) and the dilaton <f>(t) for \i = 2f3 2 = 1 

and k 2 = 0.1 . The evolution from the de Sitter bounce into a matter dominated FRW 

universe can clearly be seen. 

Our spacetime has finite curvature everywhere for t > . It can be made geodesically 

complete by patching on the contracting part of the de Sitter bounce for negative times. 

This is possible since the matching conditions for the metric are satisfied at t = 0, namely 

a(t) and a(t) are continuous. Also, <f)(t) is finite at t = and does not cause any problems. 

Therefore we have indeed obtained a nonsingular two-dimensional cosmological model. 



5 Conclusion 



Since quantum gravity is non-renormalizable[14] and at present we can extract little infor- 
mation from string theory concerning the quantum nature of spacetime, we have used a toy 
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model in 1+1 dimensions to investigate what conclusions may be reached at the classical 
level concerning singularities. The key assumption in our work is that the quantum pro- 
cesses responsible for the removal of singularities may be represented in an effective, classical 
theory of gravity, valid at curvature scales intermediate between those of any future theory 
of everything and GR. 

Motivated by the low energy limit of string theory and attempts to apply these ideas 
in four spacetime dimensions Q we have constructed a class of two dimensional models for 
gravity for which all homogeneous and isotropic solutions are nonsingular and approach 
a FRW matter dominated universe at late times. We have also demonstrated that our 
construction is valid for both the vacuum theory and, for a large class of initial conditions, 
when a conformally coupled scalar field is included in the model. 

Our construction is based on explicitly limiting the Ricci scalar (the only curvature 
invariant in two dimensions) by means of a non-dynamical scalar field <j> which we refer to as 
the dilaton. The particular choice of potential for the dilaton defines the theory and ensures 
that the solutions remain singularity free as regions of high curvature are approached. 

If the type of behaviour that we reveal here is realized by nature, there may be very 
interesting implications for cosmology. The existence of a de Sitter phase arising purely 
from the gravitational sector at early times (also seen in the four dimensional case[||) may 
provide us with a natural inflationary epoch. Although at present inflation appears to be 
the only viable solution to the horizon problem, there are difficulties associated with the fine 
tuning of parameters neccessary to implement this scenario using elementary scalar fields. 
Thus, an alternative formulation such as ours may be very attractive for cosmology. 
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Figure Captions 



Figure 1: The scale factor of the universe, as a function of the dilaton coupling, for the 
vacuum theory with an expanding de Sitter solution. 

Figure 2: An example of the scale factor, as a function of cosmological time, in the vacuum 
theory with an expanding de Sitter solution. 

Figure 3: The Ricci scalar, as a function of time, for the vacuum solution with the same 
parameters as those used in Figure 2. 

Figure 4: The dilaton, as a function of time, for the vacuum solution of Figure 2. 

Figure 5: An example of the scale factor, as a function of cosmological time, when a 
conformally coupled scalar field has been included. This example shows a de Sitter bounce. 

Figure 6: The Ricci scalar, as a function of time, for the matter solution of Figure 5. 
Figure 7: The dilaton, as a function of time, for the matter solution of Figure 5. 
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